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^ ' Abstract 
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P^ ' We propose a relationship between thermodynamic phase transitions and ground- 

\^ ■ state quantum phase transitions in systems with variable Hamiltonian parameters. 

CP . It is based on a link between zeros of the canonical partition function at complex 

temperatures and exceptional points of a quantum Hamiltonian in the complex- 
extended parameter space. This approach is applied in the interacting boson model, 
where it is shown to properly distinguish the first- and second-order phase transi- 
tions. 
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1 Introduction 



Quantum phase transitions (QPT's) are now a well documented phenomenon 
in both lattice [1] and many-body systems [2,3,4,5,6,7,8,9,10,11,12,13,14,15]. A 
QPT Hamiltonian usually reads as a superposition of two incompatible terms, 

H{\) = Ho + XV = {l-\)H{{)) + \H{l) , (1) 



[Hq, y] 7^ 0, where A G [0, 1] is a dimensionless control parameter that drives 
the system between two limiting modes of motions. It can be shown that 
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the ground-state (g.s.) average {V)o = (\E'o(A)|K|^o(A)) is a nonincreasing 
function of A. In the QPT situation it evolves in such a way that either {V)o 
itself or some of its derivatives change discontinuously (for infinite system's 
size) at a certain critical value Ac- 



Of particular interest are the cases when {V)o drops to zero and, simultane- 
ously, also {V^)o = at the critical point. Then the g.s. wave function |^o(A)) 
gets fixed for all A > Ac and {V)o plays a role of an "order parameter" whose 
value (zero or nonzero) distinguishes two quantum "phases" of the model. 
This scenario is realized in various many-body models, mostly known from 
nuclear physics. Limits H{0) and H{1) are usually connected with collective 
or single-particle motions corresponding to spherical and deformed nuclear 
shapes [2,3,4,5,6,8,9,10,11,12,13], but they can also represent paired and un- 
paired fermionic phases of nuclei [4,5,7,14], or normal and superradiant modes 
of interacting atom-field systems [15]. 



Questions often arise concerning the depth to which the term "phase tran- 
sition" can be followed in the direction of standard thermodynamics. Or is 
it only a metaphor? The g.s. QPT's in the sense of Eq. (1) happen at zero 
temperature and thus have no real thermal attributes. Nevertheless, as dis- 
cussed in Refs. [10,13,14,16], counterparts of some thermodynamic terms can 
often be derived from standard quantum-mechanical expressions. A unified 
thermodynamic-like approach to characterize the QPT situation is, however, 
missing. 



In this Letter, we present a method capable to establish the thermodynamic 
analogy for quantum phase transitions on a new, general basis. The method 
relates zeros of the canonical partition function Z{T) at complex temperatures 
T [17,18,19] to exceptional, or branch points of Hamiltonian (1) in complex- 
extended A plane [5,20,21,22,23,24,25]. This link is not artificial. While neither 
zeros, nor exceptional points occur on the real T or A axes in the finite case 
(this would imply divergences of thermodynamic observables or real crossings 
of quantum levels, which are generically forbidden), it is known that their 
distribution in the complex plane determines the key features of the system in 
the physical (real) domain. Zeros and exceptional points thus play crucial roles 
also in the fundamental theory of classical and quantum phase transitions. 
Indeed, places where in the thermodynamic limit the complex zeros of Z{T) 
approach infinitely close to the real T-axis can be identified with points of 
classical phase transitions [17,18,19], while similar convergence of exceptional 
points to real A induces singular evolution of energies and wave functions, as 
observed in T = quantum phase transitions [5,24]. 



2 QPT analog of specific heat 



One obvious way to quantitatively exploit the thermodynamic analogy relies 
on connecting the g.s. energy, Eq{X), as a function of the interaction param- 
eter A, with the equilibrium value of a thermodynamic potential, Fq{T), as a 
function of temperature T (an alternative link to the inverse temperature j3 
would not result in qualitative differences). From the relation (V^)o = dE^/dX 
it follows that if the (/ — l)-th derivative of (y)o is discontinuous at Ac, then 
derivatives of the g.s. energy are discontinuous starting from d^E^/dX^ so that 
the QPT is of the /-th order [3] . The "specific heat" defined through the sec- 
ond derivative of Eq in such a transition (in analogy to standard definition 
C = -Td^Fo/dT^) 

C, = -A^ = -X^^ = 2A ^ l(^^l^l^o)P (2) 

dX'^ dX -^Q Ei — Eq 



behaves exactly as expected for a thermodynamic phase transition of the same 
order. Here, Ei{X) and |\l'j(A)) are the i-th energy and eigenvector, respectively. 

This relation can be easily verified [13] in the interacting boson model (IBM) 
[26], where both first- and second-order phase transitions between spherical 
and deformed g.s. shapes are present in the parameter space [2,3,6,8,9,10,11,12,13]. 
The model describes shapes and collective motions of atomic nuclei in terms 
of an ensemble of N interacting s and d bosons with angular momenta and 
2, respectively. We employ a simplified IBM Hamiltonian [8,9,10,11,12,13] 



H{X) = {1-X) 



Ix-lx 



N 



+ And (3) 



where Ud = d'^ -d is the d-boson number operator, and q^ = d^s + s'^d+xid^dY'^^ 
the quadrupole operator. In the N ^ oo limit, the order parameter {V)o/N 
(normalization per boson is necessary to deal with effects of varying A^) can 
be expressed in terms of the g.s. deformation parameter (3q [2]: 

lim = ; ;77- . 4 



For X 7^ 0, the value of Po drops from a nonzero value /5oc to at A = Ac(x) = 
(4 + 2x^/7)7(5 + 2x^/7), indicating a first-order deformed-spherical QPT. For 
X = 0, the value /3o oc \/Xc — X valid in the left vicinity of the critical point 
Ac(0) continuously joins with /Sq = valid above Ac; the corresponding QPT 
is of the second order. In both cases {V^)/N ^ for A > Ac and A^ ^ oo. 



The dependence of Ci on A in the phase-transitional region is shown in Fig. 1 
for various boson numbers A^. Panels (a) and (b) correspond to x = —a/7/2 
and 0, respectively. With an increasing size of the system, the specific heat in 
panel (a), normalized per boson, would form a 5- function singularity in the 
A^ ^ oo limit (first-order QPT), while Ci/N derived from panel (b) would 
develop just a discontinuity at Ac (second-order QPT). Note that using the 
asymptotic Eq[X)/N <-> Fq{T) correspondence, one can moreover show that 
the IBM shape-phase diagram fully agrees with the classical Landau theory 
of thermodynamic phase transitions [12]. 



3 QPT distribution of exceptional points 



Exceptional points of Hamiltonian (1) are points in the complex plane of pa- 
rameter A where various pairs of eigenvalues of the complex-extended Hamilto- 
nian coalesce. They are simultaneous solutions of equations det[E — H{\)] = 
and (d/dE) det[E~H{X)] = 0, which after elimination yield the following con- 
dition [20,23]: 



D{X) = n D,{X) = - ll[E,{X) - E,{X)r = 0, (5) 

k i<j 

D,{X)= l[[E,{X)-E,{X)]. (6) 

The discriminant -D(A) is a polynomial of order n{n — 1) in A (where n is the 
dimension of the Hilbert space) with real coefficients and its roots thus occur 
as n{n — l)/2 complex conjugate pairs. Except at these points, the complex 
energy E{X) obtained from the characteristic polynomial of Hamiltonian (1) 
is a single analytic function defined on n Riemann sheets. The energy labels 
in Eqs. (5) and (6) enumerate the respective Riemann sheet according to 
the ordering of energies at real A. Exceptional points are square-root branch 
points where the Riemann sheets are pairwise connected. The leading-order 
behavior on the two connected sheets close to the branch point Aq is given 
by E{X) — E{Xo) ^ a\/X — Aq (as a doubly- valued function), with a being a 
complex constant [5,21,23]. 

The relation of exceptional points to quantum phase transitions has been 
declared several times — see, e.g., Refs. [5,9,24]. Clearly, an exceptional point 
located close to the real A axis affects the local evolution of the corresponding 
pair of real energies so that the two levels undergo an avoided crossing with 
accompanying rapid changes of wave functions. A cumulation of exceptional 
points close to some real point Ac thus can give rise to massive structural 
changes of eigenstates, as observed in QPT's. Although this mechanism was 
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Fig. 1. "Specific heat" (2) in the first-order (panel a) and second-order (panel b) 
QPT of the interacting boson model [Hamiltonian (3) with (a) x = — Vi/2 and (b) 
X = 0]. The curves, in order from the lowest to the highest, correspond to A^ = 10, 
20, 40, and 80, respectively; the insets show the increase of the maximal value with 

N. 
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Fig. 2. Exceptional points (upper panel) and crossings of energy real parts (lower 
panel) in complex A plane for J = eigenstates of Hamiltonian (3) with x = — \/7/2 
and N = 20. Only the lowest 10 levels out of n = 44 are included in the lower panel. 

illustrated by several model examples [5,24], quantitative determination of 
the arrangement and density of exceptional points needed to trigger a phase- 
transitional behavior is still missing. 



Exceptional points for levels with angular momentum J = in the vicinity 
of the IBM first- and second-order phase transitions are for A^ = 20 shown 
in the upper panels of Figs. 2 and 3, respectively. The calculation was done 
using the method described in Ref.[20]. Lower panels of both figures display — 
for the lowest ten levels, ordered according to KeEi — lines where energy real 
parts of two different levels cross, i.e., satisfy ReEi{\) = ReEj{\). As follows 
from the square-root behavior of energies close to crossings, the endpoints of 
these lines coincide with starting points of the Imii^j(A) = lm.Ej{X) lines (not 
shown here), so they represent the exceptional points. 




Fig. 3. The same as in Fig. 2, but for x = 0. In this case, single exceptional points 
(diamonds) coexist with doubly and even quadruply degenerate ones (crosses and 
rectangles, respectively). 

Several exceptional points in Fig. 3 apparently lie on the real A-axis. This 
results from the underlying 0(5) dynamical symmetry of Hamiltonian (3) for 
X = 0, with the seniority quantum number f , that survives all the way across 
the transition and thus enables the levels with the same J to cross on the real 
axis [27]. Nongeneric features implied by the additional symmetry at x = also 
include double or even higher degeneracy of some exceptional points observed 
(within the available numerical precision) in Fig. 3. 

Broad distributions of exceptional points in Figs. 2 and 3 do not seem to 
indicate the QPT behavior at A ~ 0.8, but it must be stressed that these 
figures contain points belonging to various Riemann sheets, while only the 
ground-state Riemann sheet is relevant in the present case. Unfortunately, to 
attribute the points to the respective sheets through the available numerical 
data becomes problematic already for relatively small dimensions and the 



study of effects connected with increasing boson numbers seems intractable. 



4 Alternative definition of a QPT specific heat 



The above problem can be approached via the relation of the ground-state 
related exceptional points to zeros of partition function. We will use the cor- 
respondence Z{T) ^^ -Do(^) between the partition function and k = par- 
tial discriminant (6). In fact, the square D^ of any partial discriminant is a 
polynomial with n — 1 complex conjugate pairs of roots, each of them be- 
ing simultaneously assigned to one other D^,. These roots correspond to the 
branch points located on the A;-th Riemann sheet, so that zeros of Dq rep- 
resent singularities on the ground-state Riemann sheet. A counterpart of the 
thermodynamic potential Fq = —T In Z is proportional to the g.s. "potential 
energy" , U = — J2i>o In | -Ej — -Eq I ; as obtained in the static Coulomb-gas model 
for quantum energy levels [28]. The divergence of U (zero of Dq) implies an 
infinite "force" on the ground state also in the dynamical Pechukas- Yukawa 
model [29]. For the "specific heat" one obtains (in analogy with the standard 
formula C = T^d^ In Z/dT^ + 2Td In Z/dT) 



C2 = sJ2 

i>0 



y2 I dA2 dX'2 

Ei — En 




(7) 



One can equivalently use C2 = —S\(P{\U)/d\^, which is the same as Eq. (2), 
but with Eq replaced by \U . Note that an arbitrary power of Dq used in the 
thermodynamic correspondence would modify just the overall scaling factor 
S. This factor must also depend on the size of the system (similarly as the 
thermodynamic specific heat is normalized to a mass unit) and will be dis- 
cussed later. Eq. (7) can be further decomposed [29] into a double sum of 
terms containing energy differences (in denominators) and matrix elements of 
V (in numerators), cf. Eq. (2). 

Although the quantity C2 depends solely on real-A spectral observables, its be- 
havior in a QPT is an indirect measure of the density of exceptional points on 
the ground-state Riemann sheet beyond (but close to) the real axis. Consider, 
for instance, a chain \±m = Ac ± iXm of Dq zeros along a line perpendicu- 
lar to the real axis. The -Dq polynomial is determined (up to a multiplicative 
constant) by the roots and specific heat (7) is for real A given by 

^ ,, 7pW(A'-A^),X + 2AA lJ^d\ . (8) 



where A = A — Ac and p(A) is a density of exceptional points along the \c + i\ 
line. This implies that the "latent heat" Q = \im.f.^o J^^ C2{A)dA is equal 
to zero if p(A) decreases sufficiently fast when approaching the real axis. In 
particular, if p(A) ~ A" for A -^ 0, we obtain the following possibilities: (i) a 
first-order QPT, with C2 ^ oo at A = and Q finite, for a = 0, (ii) second- 
order QPT, with C2 ^ 00 but Q = 0, for a G (0, 1], and (ni) a higher-order 
QPT, with C2 finite and Q = 0, for a > 1. This relation is the same as in 
standard thermodynamics, where the order of a phase transition refiects the 
density of the Z(T) zeros close to a critical temperature Tc [18]. It should be 
stressed, however, that in the QPT case the order deduced from C2 cannot be 
a priori expected to coincide with the order determined via Ci. We will focus 
on this problem in the following, using the IBM. 

The specific heat C2 in the IBM first- and second-order QPT is for S* = 1 
shown in Fig. 4, where panels (a) and (b) again correspond to x = —a/7/2 
and 0, respectively. All the levels with J = were included in the sum (7). We 
know that for x = some pairs of levels actually cross at real A (due to the 
seniority quantum number v) . This implies discontinuities and singularities of 
the first and second derivatives in Eq. (7), which however cancel out exactly 
and do not affect the C2 shape in Fig. 4(b). ^ 

It is clear that the peaks in panels (a) of both Figs. 1 and 4 are sharper 
and higher than those in panels (b), as indeed expected in the first- and 
second-order phase transitions. The log-log insets in these figures indicate 
that maximal values — Ai and A2 — in the Ci and C2 peaks exhibit roughly an 
algebraic increase with the boson number. The increase is faster for the first- 
order QPT than for the second-order one. We will show that with a proper 
normalization to the system's size, the values of Ai and A2 corresponding to 
the second-order QPT have a finite N —^ 00 asymptotics, while the values of 
the first-order QPT, normalized in the same way, diverge. 

Fig. 5 shows the dependence of Ai/N and A2/N'^ for the second-order QPT on 
the boson number up to A^ = 1000. The calculation for such high dimensions 
was enabled by the underlying 0(5) symmetry for x = 0- Clearly, the curves 
in Fig. 5 have finite asymptotics. The behavior of Ai is consistent with the 
analytic result Ai/N -^ 12.5 that can be derived from the N ^ 00 limit of the 
g.s. energy per boson. The faster increase of Ai in the first-order transition, 
see Fig. 1(a), indicates a divergence of Ci/N at Ac(x) for X 7^ 0- 



^ If Ei crosses with -Ej+i, the discontinuity (singularity) of the first (second) deriva- 
tives in Eq. (7) for both levels is the same, but with opposite signs, so the sum of 
both contributions is incremented correctly, as if the levels continuously passed the 
crossing. 
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Fig. 4. The same as in Fig. 1 but for "specific heat" (7) including all J = states. 

Specific heat C2 in Fig. 4 behaves in a similar way, but — as follows from 
Fig. 5 — in this case the correct normalization is by 1/N'^. This factor reflects 
the dimension of the subspace of states with J = which grows roughly as 
n ~ N'^/12 for very large boson numbers (while the total number of IBM states 
is ~ A^^/120). Since there are n— 1 pairs of exceptional points on each Riemann 
sheet, the choice of S" = l/(n — 1) oc 1/A^^ in Eq. (7) normalizes the density 
p(A) [cf. Eq. (8)] in the finite-n case to a unit integral and plays a similar role 
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Fig. 5. Normalized maximal values Ai/N (dots), 52^2/^^ (rectangles), and S2A2/N 
(triangles) of "specific heat" (2) and (7) in the IBM second-order QPT (x = 0) for 
very large boson numbers. A2 and A2 were obtained using the whole J = space or 
only its f = subspace, respectively. The behavior of A2 up to A^ = 5000 is shown 



in the inset. Scaling factors S2 
within the same range. 



10 and S2 = 0.5 were used just to show all data 



as the 1/A^ scaling of Ci in the previous case. With this normalization, the 
specific heat C2 correctly distinguishes the IBM first- and second-order phase 
transitions. 

This conclusion is also supported by the analysis of the IBM classical limit. 
Indeed, due to the degeneracy oi Pq = and Po ¥" ^ minima of the N ^ oo g.s. 
energy per boson at Ac [2], there is an exceptional point on the ground-state 
Riemann sheet that is for the infinite boson number located on the real axis, 
at the first-order QPT. Consequently, energy denominators in Eq. (7) diverge 
at this point in a way that appears to be faster than oc A^^. In contrast, there 
is no actual ground-state involving degeneracy in the asymptotic limit of the 
second-order transition, although, as can be shown, exceptional points come 
infinitely close to the real axis at the critical point. This results in a slower 
increase, which is compensated by the l/N"^ scaling. 

For X = 0, we can additionally check the validity of the S = l/{n — 1) 
normalization of specific heat (7) by considering only the v = subspace of 
J = states, with the dimension n ~ N/2. Since this subspace contains the 
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ground state and does not mix with f 7^ subspaces, one can restrict the sum 
in Eq. (7) only to these states and obtain a modified specific heat C2. It is 
expected to have a similar behavior as C2, but with S oc 1/A^. Unfortunately, 
the convergence of the maximal value A2 to some asymptotics is very slow. 
The available data for boson numbers up to 5000 (see the inset in Fig. 5) can 
only yield an upper constraint for the power in S* oc 1/N'^, namely k, < 1.35, 
which is consistent with the expectation but not fully conclusive. Work in this 
direction is in progress. 



5 Conclusions 



We have developed a method capable to indirectly measure the distribution 
of exceptional points of Hamiltonian (1) on the ground-state Riemann sheet 
near the real-A axis. This distribution plays major role in structural changes 
of the g.s. wave function. It is crucial in the determination of zero-temperature 
quantum phase transitions, but numerically difficult to be calculated directly 
for large dimensions. The method is based on the analogy between generic 
QPT arrangements of exceptional points and similar behaviors of complex 
zeros of the canonical partition function in thermodynamic phase transitions. 
The resulting formula (7) for the "specific heat" C2 normalized with respect 
to the relevant dimension, S ~ 1/n, turned out to be the appropriate measure 
for the cumulation of exceptional points near the QPT critical point. Our 
approach was tested in the interacting boson model, where C2 in the first- 
and second-order QPT's was shown to behave in the same way as standard 
specific heat in typical thermodynamic phase transitions of the respective 
orders. We expect that the method is applicable also in other QPT systems, 
even beyond nuclear physics. It discloses surprising analogy between standard 
thermodynamics and quantum mechanics of parameter-dependent systems. 
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